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Abstract 

I discuss the zeromode spectrum of lattice chiral fermions in the domain wall model 
suggested recently. In particular I give the critical momenta where the fermions cease 
to be chiral and show that the spectrum is directly related to the behaviour of the 
Chern-Simons current on the lattice. First results for the domain wall model on the 
finite lattice indicate that the relevant features of the model in the infinite system 
survive for the finite lattice. 



1 Talk presented at the International Conference on Lattice Field Theory, LATTICE 92, (1992) Amsterdam 



1 Introduction 



In this talk I want to report some first results for the domain wall model suggested by 
Kaplan [I] . The basic idea of this model is to start with an odd dimensional -and therefore 
vectorlike- theory and to add a mass term for the fermions which depends on the extra 
dimension and has the form of a soliton or kink, generating in this way a domain wall. It 
is well known that in a situation like this one finds zeromodes bound to the domain wall 
[||. Kaplan was able to show that these zeromodes represent chiral fermions on the lattice 
and that the unwanted doubler modes can be removed by introducing the usual Wilson- 
term. These results [|, which hold for the free theory and for the infinite lattice, could 
be demonstrated to survive even on a finite lattice 0. In addition the cancellation of the 
zeromode anomaly by the divergence of the Goldstone-Wilczek current as described in U 
could also been seen on the finite lattice [|. 

In this talk I will discuss the peculiar behaviour of the zeromode spectrum. I will give the 
critical momenta of these zeromodes where they lack to be chiral and relate the zeromode 
spectrum to the behaviour of the Chern-Simons current on the lattice. 



2 The infinite lattice 



To be specific I will discuss the spectrum and also the anomaly for the case of a 3-dimensional 
model. The results generalized to arbitrary dimensions can be found in ^ |5|. I start with 
the Dirac- Wilson operator on an infinite lattice with lattice spacing a = 1 

3 3 

K 3D = £ + me(s) + IJ2 A » (!) 

H=l Z n=l 

where d denotes the lattice derivative <9 M = | [5 z ,z+^ — <5z,z-^] , A the lattice Laplacian A M = 
[$z,z+n + &z,z~ii — 25 ZtZ ] , the u ^ are the usual Pauli matrices and r the Wilson coupling. I 
will denote by s the extra dimension along which the mass defect appears, while x, t are the 
2-dimensional coordinates. The domain wall is taken to be a step function e, 

r -l s < o 

e(s) = I s = (2) 
[ +1 s > 

where the height of the domain wall is given by the mass parameter m which I will choose 
to be positive throughout this paper. 

We are looking for solutions which are plane waves in the (x, t)-plane 
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where u are the eigenspinors of 03, = ±u^. With this ansatz the Dirac operator 

becomes 

2 2 r 

k sd = icr i sin (Pi) + a ^s + me(s) + r ^(cosfe) - 1) + -A s . (4) 

i=l i=l 1 

The final goal is to diagonalize the 3 dimensional Dirac operator in such a way that it reduces 
to the 2 dimensional Dirac operator for free massless fermions,^^)^/ = K 2 d^ where K 2 d 
acting on \1/ is given by 

2 

K 2 d = Yj a A = sin G°t) + °i sinfe)) . (5) 
Hence the equation to solve is 



T 

a 3 d s + me(s) — rF + -A 4 



$m ± = (6) 



where F = Y^i=t cc(l ~ cos(p$)). Following [|IJ I choose an exponential ansatz for $ away from 
the domain wall $(s + 1) = z§(s). Inserting this into @ one finds four solutions 



r - m eff ± yj m eff (m eff - 2r) + 1 

z = ! 7 

r ± 1 V ; 

where m e // = me(s) — rF, the ± in the nominator stand for the two roots and the ± in the 
denominator stand for the chirality. Note that in the limit r = 1 eq.(^) can be reduced to 
the corresponding expressions in [|l[ . 

One has to impose the condition that the solutions are normalizable to obtain sensible 
wavefunctions. This means that \z\ > 1 for s < and \z\ < 1 for s > 0. The boundaries of 
the regions where chiral solutions exist are obtained by setting \z\ = 1. Explicit matching 
of the normalizable solutions for positive and negative s at s = enables one to determine 
the regions with chiral fermions. One finds that existence and chirality of the solutions is 
independent of the sign of r and that a negative m leads to opposite chiralities. Depending 
on the values of m/r one gets m = rF and m = r(F + 2) as the boundaries for the critical 
momenta, where F is defined as above. 

The results can be summarized as follows. Starting with m/r = one finds no chiral 
fermions. For increasing < m/r < 2 the region in momentum space around p = (0,0) 
where chiral modes exist grows. This region is bounded by m = rF which gives the upper 
critical momenta. Increasing m/r above m/r = 2 opens up the two "doubler" modes at 
p = (0, 7r) and p = (ir, 0) which have flipped chirality, while the original mode at p = (0, 0) 
disappears. Here the boundaries of the regions in momentum space are given by m = rF for 
the lower and m = r(F + 2) for the upper critical momenta. For m/r > 4 the two "doublers" 
disappear and one gets a zero mode at p = (n, it) with the same chirality as the mode at 
p = (0,0). The boundary for the lower critical momenta is given by m = r(F + 2). This 
mode is finally also lost as m/r is increased to m/r > 6. 

It should be remarked that this spectrum stems from solutions only, and that there 
are no solutions for positive m. The change of the chirality is the usual reinterpretation 
of the chirality at different corners of the Brillouin zone. 
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The zero mode spectrum as found here is directly related to the coefficient of the lattice 
Chern-Simons current induced by heavy fermions. This current is responsible for the anomaly 
cancellation. As there is a chiral zeromode bound to the domain wall, there exists the 
corresponding anomaly. However, we started with a vectorlike theory and consequently the 
theory should be anomaly free. As is explained in [|l|] this contradiction is resolved by the 
fact that the domain wall induces a Goldstone-Wilczek current |J far off the domain wall 
the divergence of which exactly cancels the zeromode anomaly on the domain wall @. 

Recently, this current was also calculated on the lattice. The radiatively induced Chern- 
Simons action in d=3 dimensions is given by || 



cs 



d 3 xA„d„A p . (8) 



The effective action is then given by cTcs- The coefficient c can be calculated in perturbation 
theory || || and is given by 

C = itfiup I\ q =o (9) 



f d 3 p 

I ?2^)3 Tr]iS ^ A ^ p,p ~ q " > 5 (P~<7)Ap(p + <?,p)] (10) 



d{q) u 

with / the following integral 

d 3 p 
(2^ 

where S is the fermion propagator and A M denotes the photon vertex. Imposing the Ward 
identity 

A^p,p) = -id/d Pfl S- 1 (p) (11) 

the relevant integral becomes 

d 3 p 



I ttt^Tt \Sd.S- 1 } \SdvS- 1 } \SdpS- 1 } . (12) 



(2tt)* 

The fermion propagator can be written in a generic form, S^ 1 (p) = a(p)+ib(p)a which has 
the structure of S~ 1 (p) = N(p)V(p), where V(p) is a SU(2) matrix and N(p) a normalization 
factor. Provided that S~ l does not vanish for any p the integral eq.(0) does not depend on 
N(p) and has a simple topological interpretation: It is the winding number of the map V 
from the torus T 3 to SU(2) or the sphere S 3 . (In general this map is a map from the torus 
T d to the sphere S d .) 

Taking the usual Wilson fermion propagator the winding number can be calculated by 
using the fact that the integral has to be computed only in infinitesimal regions near the 
Brillouin corners. One obtains (see || for more details) 

c -^h-^) k ^^h- as) 

327T 5 ^ \m — 2rk\ 

This expression should be compared with the (unregulated) continuum result 

—im 

c cont 777^ si i" • l-L^V 

627r 6 to 
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Therefore one finds that for m/r < 0, c = 0. For m/r < 2 the lattice result is twice the 
continuum value. For 2 < m/r < 4 it is — 4c con4 and for 4 < m/r < 6 it is again twice the 
continuum value. 

For the case of the domain wall model in which we are interested here the result has the 
following implication: For the calculation of the Goldstone-Wilczek current one goes far off 
the domain wall and assumes there the mass to be constant. Then one can calculate the 
current with the method described above. For m/r < 2 one finds therefore that the current 
only flows on one side of the domain wall and has twice the value of the continuum result. Of 
course, the divergencies will come out the same. The behaviour of the lattice Chern-Simons 
coefficient finds its exact correspondence in the zeromode spectrum as discussed above where 
we found 1 lefthanded, 2 righthanded and 1 lefthanded chiral fermion for < m/r < 2, 
2 < m/r < 4 and 4 < m/r < 6, respectively. 



3 And the Finite Lattice 



Any numerical work on this system will necessarily involve finite lattices, and so I now 
compare the results obtained on the infinite system with the ones of a finite lattice. On 
the finite lattice one has to choose some boundary conditions. Taking periodic boundary 
conditions generates a second anti-domain wall. The mass term is therefore modified to be 
mei(s) with 

r -l 2 < s < if 

e L (s) = \ +1 ^ + 2 < s < L s . (15) 
I s = 1,^ + 1 

The zeromodes on the finite lattice can be searched for by solving the Hamiltonian 
problem numerically. If one again assumes plane waves in the x-direction the Hamiltonian 
is given by H 



r 

H = -at [icr 2 sm[p x ) + cr 3 d s + me L (s) + r(cos(p x ) - 1) + -A s 



(16) 



The eigenvalues and eigenfunctions of the Hamiltonian eq.(|16l) were calculated numeri- 
cally . To find the critical momenta the ratio 

^t— r ( 17 ) 

was studied, which is a normalized measure for whether the fermions are chiral or not. It is 
zero if the fermions are chiral and R > for non-chiral modes (see fig. 2b in ||). To determine 
whether one still has chiral fermions a threshold value for R was defined. If R < 0.01 the 
fermions were regarded to be chiral. 

Comparing the results from the infinite system with the finite lattice calculations with 
L = 100 P], one finds that the two are practically indistinguishable. For L = 20, a lattice 
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size realistic for simulations, a small shift occurs. Fixing m and r we find for m/r < 2 a 
smaller value and for 2 < m/r < 4 a larger value of the critical momentum. 

It is also possible to extract the Chern-Simons current on the finite lattice in the presence 
of a smooth external gauge field configuration. The current is most easily computed by 
the inverse fermion matrix of the model which can be obtained by standard methods like 
conjugate gradient. Note that this is an exact solution of the numerical problem and that 
no simulation is involved 01] . 

I show in fig. 1 the Chern-Simons current on a 16 3 lattice with m = 0.81 and r = 0.9. 
Note that these values of m, r are quite large. The picture nicely demonstrates that the 
advocated flow of the current as obtained above is reproduced on the finite lattice: It flows 
only on one side of the domain wall and is zero on the other. 

It is a simple task to get the divergence of the current and one can demonstrate @ that 
it obeys the continuum anomaly equation to a very good accuracy. It is also possible to 
find non-trivial anomaly cancellation like in the 3-4-5 model where the individual fermion 
currents are anomalous but the sum of them cancel. 

In summary the domain wall model shows a lot of promising features on the finite lattice 
like the chiral zeromode spectrum and the correct anomaly behaviour. Although these prop- 
erties were only found for the free theory or with weak external gauge fields they certainly 
point into the right direction and are encouraging. The next step is to include dynamical 
gauge fields and see whether one can find the desired behaviour. 
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Figure Caption 

Fig.l I show the Chern-Simons current in arbitrary units as obtained on a 16 3 lattice. It 
shows the expected peculiar behaviour that it flows only on one side of the domain wall and 
is zero on the other. The locations of the domain walls is at s=l and at s=9. 
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